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Universal Fermi gases in mixed dimensions 
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We investigate a two-species Fermi gas in which one species is confined in a two-dimensional plane 
(2D) or one-dimensional line (ID) while the other is free in the three-dimensional space (3D). We 
discuss the realization of such a system with the interspecies interaction tuned to resonance. When 
the mass ratio is in the range 0.0351 < m2n/m^, < 6.35 for the 2D-3D mixture or 0.00646 < 
miD/m3D < 2.06 for the 1D-3D mixture, the resulting system is stable against the Efimov effect 
and has universal properties. We calculate key quantities in the many-body phase diagram. Other 
possible scale-invariant systems with short-range few-body interactions are also elucidated. 



Introduction. — Experiments using ultracold atomic 
gases have attracted considerable interest because of 
their high designability and tunability. Not only can 
the strength of interaction be varied via the Feshbach 
resonance, but also the dimensionality of space can be 
changed by means of strong optical lattices. The Bose- 
Einstein condensation in trapped atoms has been realized 
both in one and two dimensions A strongly inter- 
acting one-dimensional Fermi gas has been observed Q, 
which is an atomic realization of Tomonaga-Luttinger liq- 
uid. 

What the cold atom community has not paid atten- 
tion to is the system in mixed dimensions, namely, the 
system where different types of particles live in different 
spatial dimensions. Actually the idea of mixed dimen- 
sions is commonly used in various subfields of physics. 
For example, in the brane world model of the Universe, 
ordinary matter is confined in a three-dimensional space 
embedded in higher dimensions where gravity can propa- 
gate 0] • Also, recently realized graphene can be regarded 
as a system in mixed dimensions where photons induce 
the three-dimensional Coulomb interaction between elec- 
trons confined in a two-dimensional sheet 

In this Letter, we discuss the realization of an analo- 
gous system using cold atoms with two fermionic species 
A and B, where A atoms are confined in a two- 
dimensional plane (2D) or one-dimensional line (ID) 
while B atoms are free in the three-dimensional space 
(3D). At low energies, we show that the interspecies in- 
teraction in a certain range of mass ratio is characterized 
by a single parameter, the effective scattering length a e g, 
whose value is arbitrarily tunable by current experiments. 
In particular, if the effective interaction is tuned to reso- 
nance |a c ff| = oo, such a system is scale invariant and has 
universal properties. (Here intra-species interactions are 
neglected and we use the term "universal" in the narrow 
sense that the interaction does not have any dimension- 
ful parameters.) The corresponding many-body system 
forms a novel type of unitary Fermi gas in mixed dimen- 
sions. 

Experimental realization. — To confine only A atoms, 
we note that the strength of the optical trap experienced 
by neutral atoms is proportional to the polarizability 



which depends on the atomic species and the laser fre- 
quency lot,. By tuning u>t, at a zero of the polarizability 
of B atoms, which exists between any two nearest res- 
onances, we can confine only A atoms. If the optical 
trap with such a laser frequency is d± dimensional with 
d± = 1 or 2, A atoms are confined in 2D or ID, respec- 
tively, while B atoms are free from the confinement. 

If the confinement potential applied to A atoms is har- 
monic with the oscillator frequency lj±_ , an A atom inter- 
acting with a B atom is described by the Hamiltonian: 

H = ~^kr + 5 m A"x x AX - + V( x a,x b ). We 

introduce a notation xa = (xa\\ , xa±), where xa± is the 
d± -dimensional coordinates affected by the confinement 
potential. V(xa, Xb) is the bare interspecies interaction. 
In the limit of zero-range interaction, it is characterized 
by a single parameter, the s-wave scattering length a, 
which is arbitrarily tunable by means of the Feshbach 
resonance. 

Because of the partial confinement of the A atom, the 
scattering property between A and B atoms is modified 
from the free-space scattering. In particular, an infinite 
number of confinement-induced resonances appears. To 
determine the position of resonances as a function of l±_ /a 

with l± = m ^ ± , we solve the Schrodinger equation 

H'^(xa 1 xb) — E^{xa,xb) in the zero-energy limit E — 
4rftw_L — > 0. After excluding the center-of-mass motion 
in the parallel direction, the solution is written as 

* = J d d± y ± G(x A \\ - x B \i,XA±,x B x;0,yx,y±)f(y±) 

+ C<t>a{x A x), (1) 

where 4 > oi x A±) — ex p[ — x a±/(^D] * s * ne g roun d state 
wave function of the d^-dimensional harmonic oscillator. 
G is the Green's function given by 
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with u = rriA/rnB being the mass ratio. The un- 
known function f(y±)/C is uniquely determined by im- 
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FIG. 1: (left) a e g/l± as a function of l±/a for d± = 1 (2D- 
3D mixture). The mass ratio is chosen to be u — 0.15. The 
vertical dotted line indicates the position of the broadest res- 
onance, (right) Lowest five values of Zj_/a to achieve the res- 
onances |a e ff|/Zj_ = oo as functions of u. 



posing the short-range boundary condition ^(xa^xb) oc 
\* a ^* b \ -Z + °(\ x a- x b\) when \xa-x b \ -> [|. The 
asymptotic form at a large separation \xa — Xb\ — ► oo, 
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^o(^a±), (3) 



gives an effective scattering length a c ff, and whenever 
a c ff diverges, a resonance occurs. Here T>(xa,x b ) = 

with TTlAB 



being 



niAB B± ""/1J3 vriA+niB 

the reduced mass is the "interparticle distance" with the 
anisotropic weights to the parallel and perpendicular di- 
rections. 

Numerically obtained a c g/l± as a function of l±/a is 
plotted for d± = 1 (2D-3D mixture) in Fig. Q] and for 
e?^ = 2 (1D-3D mixture) in Fig. [2j We chose the mass ra- 
tio u = 0.15 corresponding to the physical case of A = 6 Li 
and B = 40 K. These figures show that a e ff can be tuned 
to any desired value by varying a or l±. In particular, the 
values of l±/a to achieve the resonances |a c ff|/Zj_ = oo are 
plotted as functions of u in Figs. Q] and [2] 

If the confinement length l± is much smaller than any 
other length scales of the system such as a e g and a mean 
interparticle distance at finite density, one can neglect 
the motion of A atoms in the confinement direction. 
Consequently, the resulting system is mixed-dimensional, 
namely, A atoms are confined in 2D or ID while B atoms 
are in 3D. The interspecies interaction is solely charac- 
terized by the effective scattering length a c g defined in 
Eq. When a e g > 0, there is a shallow bound state 
with the binding energy ^binding = ~ ^/(^ABOeff)- 

If the effective interaction is tuned to resonance |a e ff| = 
oo, the system is scale invariant and universal, i.e., inde- 
pendent of the short-range physics. The corresponding 
many-body system forms a novel type of unitary Fermi 
gas in mixed dimensions. In the rest of this Letter, we 
concentrate on the system at the resonance. 

Stability of the system at resonance. — An important 
question regarding our novel system is its stability. In 
the usual 3D case, it is known that two heavy and one 
light fermions with mass ratio greater than 13.6 develop 
deep bound states in the limit of zero-range interaction 



FIG. 2: Same as Fig.[T]for d± = 2 (1D-3D mixture). 

(Efimov effect) jfjj]. Accordingly the corresponding many- 
body system is not stable toward collapse. To establish 
the stability of the unitary Fermi gas in mixed dimen- 
sions, we study three-body problems and determine the 
range of u where the Efimov effect is absent. 

We first consider the problem of two A atoms in- 
teracting with one B atom. The wave function 
$?(xai, XA2t Xb) for V(xAi,XB) > satisfies the 



* = 



Schrodinger equation: 

° ^ 2mA 2mA 2mi 

. Now xa = Xa\\ is two- or one-dimensional coordi- 
nates. The resonant interspecies interaction is taken into 
account by the short-range boundary condition: ^> oc 
V{x A i,x B Y l + 0[D(xAi,x B )] when V(x A i,x B ) -> 
[see Eq. In hyperspherical coordinates, the wave 

function at a short distance behaves as 0] 



i£(xai,xa2,xb) ^ R 7 tpi(fl), 



(4) 



where R is the hyperradius and Q denotes the hyperan- 
gular variables. 7 is the scaling exponent of the wave 
function classified by the quantum number / and can be 
determined by writing ^ in terms of the Green's function 
and imposing the short-range boundary condition [8] . 

In the 2D-3D mixture, I — 0,±1,±2, ... is the orbital 
angular momentum projected to the 2D plane and 7 sat- 
isfies the equation [9|] 
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with 



3(7, 0) = ^ r(- 7 - 1) T( 7 + |) e-^+V 
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Here 2F1 is the hypergeometric function. Equation ([5]) 
has a pair of solutions, 7+ and 7_, related to each other 
by 7+ + 7_ = —3. For the p-wave channel \l\ = 1, as the 
mass ratio u increases from to the critical value Umax = 
6.35, 7+ decreases from to — |. When u is increased 
further, u > u max , j± become complex conjugates y± = 
— I ±zs ; indicating that the wave function ([J| oscillates 
as Rr 3 / 2 sin(so In R+5) and the Efimov effect takes place. 
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The critical mass ratio for other odd I is greater than 6.35 
and the Efimov effect is absent for any even I. 

On the other hand, in the 1D-3D mixture, / = ±1 in 
Eq. (0| corresponds to the parity of the wave function 
and 7 satisfies the equation [1 01 ] 



y/2u + 1 _ cos [(7 + I) a] + Zcos[(7 + 1) (tt - a)] 
u + 1 (7 + 1) sin [tt (7 + 1)] 



(7) 



with a = arccos ■ Similarly, this equation has a 
pair of solutions, 7+ and 7_, related to each other by 
7 + +7- = —2. For the odd-parity channel I = — 1, when 
u exceeds the critical value u max = 2.06, 7± become com- 
plex conjugates j± = — 1 ± iso, indicating the Efimov 
effect. In the even-parity channel, the Efimov effect does 
not take place. 

We then consider the problem of one A atom in- 
teracting with two B atoms. The wave function 
^(xa, xbi, XB2) — TOipiip) satisfies the Schrddinger 
equation and the short-range boundary condition anal- 
ogous to the previous case. 7 in this case satisfies an 
integral equation Q ■ For the p-wave or odd-parity chan- 
nel, we obtain 
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(8) 

where p = (p\\,p±.) and S p = S (|p_L.|/|p|i|) is an un- 
known function. There is a pair of solutions, 7+ and 7_, 
related to each other by 7+ + 7- = —4. As the mass 
ratio u decreases from 00 to the critical value u m i n , 7+ 
decreases from to —2. When u is decreased further, 
u < w mim 7± become complex conjugates y± = — 2±iso, 
indicating the Efimov effect. The critical mass ratios are 
found to be w m in = 0.0351 for the 2D-3D mixture and 
Wmin = 0.00646 for the 1D-3D mixture. In other chan- 
nels, the Efimov effect does not take place for the mass 
ratio greater than u m [ n . 

So far we have supposed that both atomic species A 
and B arc fcrmionic. If either A or B atoms are bosonic, 
one can confirm that the Efimov effect takes place for 
any mass ratio in the s-wave or even-parity channel. 
Therefore, for the stability of many-body systems, both 
atomic species have to be fermionic with the mass ratio 
between M m ; n and w max . For example, the combination 
of fermionic atoms, A = 6 Li and B = 40 K [u = 0.15), 
can be used to realize both 2D-3D and 1D-3D mixtures, 
while the opposite combination, A = 40 K and B = 6 Li 
(u = 6.67), suffers the Efimov effect in both mixtures. 

Many-body physics. — We now proceed to the many- 
body physics of unitary Fermi gas in mixed dimensions. 
It is convenient to work in the grand-canonical ensemble 
by introducing chemical potentials \ia to A atoms and 
[iB to B atoms. The phase diagram of the system in 
the plane of \xa and /ib is depicted in Fig. [3] There are 
four distinct regions: vacuum (/i^ < 0, [is < 0), pure 
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FIG. 3: Phase diagram in the plane of /j,a and fiB at zero 
temperature. 



2D or ID gas {^a > 0, [ib < Ms'); P ure 3D gas (fiB > 
0, [i a < MX); and mixed gas (rest of the phase diagram). 
As we discuss below, the mixed-gas phase can be divided 
into two phases: mixed gas I (/zb < 0) and mixed gas II 
(/is > 0). It is possible that these mixed-gas phases are 
further divided into several phases by nontrivial many- 
body physics including inter- and intra-species pairings. 

Because of the lack of translational symmetry in the 
perpendicular direction, the density of B atoms ns(|a;^|) 
is a nontrivial function of the distance from the 2D plane 
or ID line. If /is < (mixed gas I), all B atoms are 
localized around the 2D plane or ID line because they 
are attracted by A atoms. Accordingly B atoms have a 
vanishing density at infinity ns(|a;j_| — > 00) — > 0. When 
Hb exceeds zero (mixed gas II), a portion of B atoms 
are no longer bound due to the Fermi pressure. In this 
phase, the density of B atoms at infinity is given by that 

x (2 mB ^ B f 2 



of noninteracting fcrmions Ub(\xa 



In order to make the discussion more quantitative, we 
estimate the key parameters in the phase diagram. The 
chemical potential at the phase boundary H°a\b) normal- 
ized by the Fermi energy of majority atoms Hb(a) (cor- 
responding to the slope in Fig. [3]) is a universal number 
because of the scale invariance of the interaction. To cal- 
culate /i^'Y/xs, we employ a many-body variational wave 
function assuming single particle-hole excitations 111 ]: 



\q\>hk FB 

\X) =Xo|0||)jFS) B + Xp, q \P\\-q\\) A \p,q) B - (9) 

\p\<hk FB 



Here \p\\) a is the momentum eigenstate of a single A 
atom and |FS) B is the unperturbed Fermi sea of B atoms 
with the Fermi momentum hkps = ^ZmsfJ-B- \PiQ)b 
describes the perturbed Fermi sea with the particle (mo- 
mentum q) and hole (p) excitations. This simple ap- 
proximation is known to give reasonable results in the 
usual 3D case 12]. By minimizing the energy expecta- 
tion value with respect to the coefficients xo and Xp,qi 



4 



n.o r 



0.5 1.0 2.0 5.0 



oni om 005 n io o?n 50 l.oo 2.00 



-0.5 
-1.0 
-1.5 

-2.0 / ti"'/ti B 
-2.5 
-3.0' 




FIG. 4: I ' \ib and u^/fiA as functions of u for the 2D-3D 
mixture (left) and 1D-3D mixture (right). 



we obtain the equation to determine /i™: 
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(10) 



Similarly, it is straightforward to derive the equation to 
determine fig 1 . Numerically obtained critical chemical 
potentials fi c j^ / fig and fi c g / fiA a s functions of the mass 
ratio w m in < u < u ma x are plotted in Fig.[l]for the 2D-3D 
and 1D-3D mixtures. 

Other universal systems. — Finally we elucidate other 
possible scale-invariant systems with short-range few- 
body interactions. First it is possible for an additional 
intra-species interaction to produce a resonance in a 
three-body channel [13j . AAB three-body resonance can 
be introduced if the three-body wave function (J4| with 
7 = 7_ is normalizable at origin R — * 0. This corre- 
sponds to the range of mass ratio 2.33 < u < 6.35 for the 
2D-3D mixture or < u < 2.06 for the 1D-3D mixture. 
Similarly, ABB three-body resonance can be introduced 
if the mass ratio is in the range 0.0351 < u < 0.0661 
for the 2D-3D mixture or 0.00646 < u < 0.0202 for the 
1D-3D mixture. The corresponding many-body systems 
become universal Fermi gases with both two-body and 
three-body resonances. 

It is also possible to realize scale-invariant systems in 
different spatial configurations from those considered so 
far. The key idea is to find a system in which the relative 
motion is characterized by three coordinates (excluding 
the center-of-mass translations if exist). In such a system, 
the interspecies short-range interaction is solely charac- 
terized by the effective scattering length and at resonance 
the scale invariance is achieved. Consider a 1D-2D (2D- 
2D) mixture where B atoms are confined in a 2D plane 
having a point-like (linear) intersection with the ID line 
(other 2D plane) in which A atoms are confined. Here 
A-B relative motions are three-dimensional and the re- 
sulting systems at resonance provide two more classes of 
universal two-species Fermi gases in mixed dimensions. 

Our idea can be extended to cases with more atomic 
species. If the number of atomic species is three (A, B, 
and C), there exist two classes of universal three-species 
Fermi gases in mixed dimensions. One is a 1D 2 -2D mix- 
ture where A and B atoms are confined in the same ID 



line embedded in a 2D plane where C atoms are confined. 
The other is a 1D-1D-1D mixture where A, B, and C 
atoms are independently confined in ID lines and these 
three lines intersect at a single point. In cither case, the 
interspecies three-body interaction is three-dimensional 
and can be in principle tuned to resonance. 

Furthermore, if we have four atomic species and all of 
them are confined in the same ID line, the interspecies 
four-body interaction becomes three-dimensional. The 
corresponding many-body system at resonance will be 
especially interesting because it forms a universal four- 
species Fermi gas purely in one dimension. 

Conclusion. — In addition to the well-known s-wave 
two-body resonant interaction in pure 3D, there exist 
seven more types of scale- invariant short-range few-body 
interactions in pure and mixed dimensions: two-body 
resonant interactions in 2D-3D, 1D-3D, 2D-2D, and 1D- 
2D mixtures; three-body resonant interactions in 1D 2 -2D 
and 1D-1D-1D mixtures; four-body resonant interaction 
in pure ID. At finite densities, each of them corresponds 
to a novel class of universal multispecies Fermi gases that 
have nontrivial many-body physics. We have given de- 
tailed analyses in the cases of two-species Fermi gases in 
the 2D-3D and 1D-3D mixtures. 
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